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ABSTRACT 

7-ray radiation from pulsars is usually thought to be mostly produced by the synchro¬ 
curvature losses of accelerated particles. Here we present a systematic study of all cur¬ 
rently reported, good-quality Fermi-hAT pulsar spectral data. We do so by applying a 
model which follows the particle dynamics and consistently computes the emission of 
synchro-curvature radiation. By fitting observational data on a case by case basis, we 
are able to obtain constraints about the parallel electric field, the typical lengthscale 
over which particles emit the bulk of the detected radiation, and the number of in¬ 
volved particles. The model copes well with data of several dozens of millisecond and 
young pulsars. By correlating the inferred model parameters with the observed timing 
properties, some trends are discovered. First, a non-negligible part of the radiation 
comes from the loss of perpendicular momentum soon after pair creation. Second, the 
electric field strongly correlates with both the inverse of the emission lengthscale and 
the magnetic field at light cylinder, thus ruling out models with high-energy photon 
production close to the surface. These correlations unify young and millisecond pul¬ 
sars under the same physical scenario, and predict that magnetars are intrinsically 
7-ray quiet via syncrhro-curvature processes, since magnetospheric particles are not 
accelerated enough to emit a detectable 7-ray flux. 


1 INTRODUCTION 

Thanks to the increasing number of detections by Fermi- 
hAT, the number of known 7 -ray pulsars is currently in ex¬ 
cess of 161 Q 71 of them are millisecond pulsars (MSPs), 
with periods between 1.5 and 22 ms, while the remaining 
90 are young pulsars (YPs), with periods up to half a sec¬ 
ond. Data for 117 of these pulsars were presented in de- 
tail in the second Fermi-hAT pulsar catalog (2PC hereafter, 
lAbdo et al.l 1201^ 1. We shall draw the targets of our study 
from this impressive database. These pulsars are among the 
most energetic, and can be quiet or loud in radio and/or 
X-rays. Systematic studies of the 7 -ray pulsars usually re¬ 
gard their pulse profile, and the constraints it can impose 
on the magnetospheric geometry under some particular as¬ 
sumptions of a given model for its location llRomanil 1 19961 : 
iRoniani fc Wattersll 201 (]l : IPierbattista et al.ll2015l ). Here, we 
present a systematic study and fit the phase-averaged spec¬ 
tra of the 81 pulsars having good-quality data in the 2PC. 

This works builds upon the result of a research pro¬ 
gram about the high-energy synchro-curvature (SC) ra¬ 
diation coming from pulsar magnetospheres. Within this 
program, we have earlier reformulated the SC expressions 
and computed particle trajectories under this kind of loss 


dVigano et al.|[^15a^ . We have also studied the assumptions 
and caveats in the semi-analytical thin and thick outer gap 
(OG) models dVigano et al.ll^l5bl l3l. Finally, we have pro¬ 
vided phase-averaged and phase-resolved analysis of the 7 - 
ray s pectra of Crab, Geminga, and Vela dVigand fc Torre^ 
I2OI5II . The results obtained in the latter work have prompted 
us to undertake this systematic analysis. Our aim here is to 
test whether the expected SG radiation is compatible with 
the 7 -ray spectra of the Fermi-LAT detected pulsar popula¬ 
tion, and if so, to obtain general conclusions and explore cor¬ 
relations between model parameters and observational data. 
We tackle these aims by fitting the phase-averaged spectra 
of a large population of YPs and MSPs. 


In ^we summarize the theoretical model used and how 
we predict the observed spectra (see also Appendix m and 
references therein), focusing on the only three parameters 
used in the spectral fits. In we define the sample and 
present the results of the spectral fitting for each of the 
pulsars (see also Appendix H). m is dedicated to the explo¬ 
ration of correlations between the inferred physical param¬ 
eters of the best-fitting models and the timing and pulse 


^ This number of pulsars appears already in the list at 
https://confluence . slac. Stanford. edu/display/GLAMCOG/Public+Li4tFtl}fl£Li5'C!f@5®4^(@<0aMfe£RHy-@iIIli4®e In ^ we discuss the 

updated November 2014. implications of these results. 
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Table 1. Parameters entering in the SC model, definition, range of observed/expected values, reference to further discussion justifying 
the adopted range or definition, and treatment in the fitting procedure. The first block of parameters have directly observed values, 
the second block lists the only three parameters varied in the fitting procedure, and the third block collects the less spectra-affecting, 
virtually unconstrainable parameters. 


Parameter 

Definition 

Range 

Ref. 

Treatment 

p 

Spin period 

1-500 ms 

11 ] 

measured 

p 

Spin period derivative 

0 

1 

to 

r 

0 

1 

to 

[ 1 ] 

measured 

E|| 

Parallel electric field 

10®-10“ V/m 

[ 2 ] 

fit 

xo/R\c 

Lengthscale of the bulk 7 -ray emission, Eq. |(2j 

0 

1 

1 

0 

0 

[ 2 ] 

fit 

No 

Effective number of particles, Eq. |[3 

1027_io34 

[ 2 ] 

fit 

V 

Radius of curvature position-dependence, rc{x) = Ric {x/R\c)^ 

0 . 2 - 1.0 

[3] 

fixed to 0.5 

b 

Magnetic field position-dependence, B{x) = Bs (Rir/x)^ 

2-3 

[3] 

fixed to 2.5 

^in/^lc 

Inner gap location 

0 . 2 - 1.0 

14] 

fixed to 0.5 

/ R\c 

Outer gap location 

1 . 0 - 2.0 

14] 

fixed to 1.5 

Fin 

Particle Lorentz factor at birth 

103-10'‘ 

[4] 

fixed to 10 ® 

ttin 

Particle pitch angle at birth 

1 

0 

[4] 

fixed to 7 r /4 


References; | 1 ] 2PC; [ 2 ] §2.3 of IVigano fc Tor^ (l2015h ; [3] Appendix of IVigano et al.l (l2015lJl : |41 IVigano et al.l ll2015cll . 


2 THE SYNCHRO-CURVATURE MODEL 

The model, presented in detail in a series of pre vious works 
dVigano et alj[^15al f9l3: IVigand fc Torresll2015l l. calculates 
the SC spectra emitted by a particle accelerated along a gap 
placed in some region of the pulsar magnetosphere (within 
or outside the light cylinder). We refer the reader to the 
m ost important formul a e, quoted in the Appen dix lAl and 
to lVigand et al.l l|2015al l: IVigano fc Torrei l|201lTh . which no¬ 
tation we follow, for further explanations, while here we shall 
only sketch the main conceptual points. 

We assume the existence of a parallel component of the 
electric field in a gap, which is responsible of accelerat¬ 
ing particles up to high Lorentz factors, F. We calculate the 
evolution of the kinematic properties of said particles while 
accelerating, F and the pitch angle a, by numerically solving 
the equations of motion along the particle movement when 
the full SC losses are taken into account. Together with the 
local (i.e., varying along the gap extent) values of the mag¬ 
netic field B and of the radius of curvature F and a 
regulate the local emission of SC radiation. The evolution 
of these parameters along the particle trajectory is consis¬ 
tently followed when computing the radiation-reaction force 
(given by the SC power) and the spectral energy distribution 
of the radiated emission, dPsc/dEj. The latter is a cumber¬ 
some expression (see Eqs. A2 to AlO of Appendix [All but 
can be numerically computed. 

Once we have computed the particles trajectories for 
each given set of parameters, we assume an effective parti¬ 
cle distribution, dNe/dx{x\xo, No,Xin,Xout), along the trav¬ 
elled distance on the magnetic field line, x, between the as¬ 
sumed boundaries Xin and Xout (which are parameters in the 
model, see below). Thus, in order to obtain the detectable 
radiation we convolve the single-particle SC spectrum with 
dNe/dx: 

dPgap dP„, dN,, , 

The effective particle distribution represents only the num¬ 
ber of particles, per unit of distance, emitting radiation to¬ 


wards us. For representing the particle distribution we shall 
assume 

J/V 

_ jv ^ _ ( 2 ) 

dx ° a;o(l — e“^°“‘Po) 

This functional form can be used to weight differently (ac¬ 
cording to the value of so) the number of particles at dif¬ 
ferent places of the gap extent. Put otherwise, it allows to 
study whether the synchrotron-like emission (i.e., the loss of 
the initial component of the particle momentum perpendic¬ 
ular to the magnetic field line, which happens close to the 
pair creation location), or the curvature-like emission (i.e., 
the loss of the parallel component of the particle momen¬ 
tum, dominating at larger distances) dominates the emitted 
spectrum. 

By expanding Eq. 0 above using the details found in 
Appendix one can see that the model contains eleven 
parameters. The model parameters are, however, neither all 
unknown, nor equally relevant. Two of the parameters are 
In fact the observable spin period P and its derivative P, 
which define in turn the light cylinder radius, defined as 

Ric = = 4.77 X lO^Pfs] cm , (3) 

27r 

and the inferred value of the magnetic field at the polar 
surface, Bs (see iMl). Then, six other parameters regulate 
the particle trajectory: E\\, 77, b, Xi^, Fin, and ain; and three 
additional parameters describe the effective particle distri¬ 
bution: No, xq/R\c, and XoMt/R\c- Table [T] gives details on 
these parameters, their considered range, and provide spe¬ 
cific re ferences to further discu ssion. 

In IVigano &; Torres! ll2015ll we proved that only three 
parameters are the most relevant in determining the spec¬ 
trum: the parallel electric field in the gap, Ey, the normaliza¬ 
tion of the effective particle distribution, Nq, and the weight¬ 
ing parameter so- The other model parameters barely affect 
the predicted spectral outcome if they vary within the con¬ 
sidered range, and we fix them here to fiducial values for all 
pulsars, as Indicated in the third block of Table [T] Those six 
fixed parameters, which will not be discussed further here, 
are the initial Lorentz factor Fin, the initial pitch angle ain, 
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Table 2. Timing parameters, timing-inferred properties, 0.1-100 GeV luminosity, and best-fitting parameters xq and A^o) to the 

phase-average spectra for the 59 YPs of our sample. We indicate lower (upper) limits for E’jj (xq and xq/Ric) in the cases with strong 
degeneracy (see text and Appendix Ell. When distance is unknown, we use its upper limit (i.e., the distance to the edge of the galaxy 
in that particular direction) to have an upper limit on L, which is used to constrain A^q (which, therefore, has to be taken as an upper 
limit). The last columns give the value of for the best fitting model and the number of considered data points (the number of degrees 
of freedom (dof) is A^bin ~ 3). 


Pulsar 

P 

log P log E 

log T 

log Bs 

log Sic 

log L 

log Bu 

log IQ 


log A^’o 

^min^-^bins) 


[m=] 

[erg/s] 

[y] 

[G] 

[G] 

[org/B] 

[V/m] 

[cm] 




J0007 + 7303 

315.9 

-12.45 

35.65 

4.15 

13.04 

3.50 

34,97 

> 

7.87 

< 

7.09 

< 

-2.09 

01 QC+0.01 

■^I'-^b-o 03 

41.37 (11) 

J0106+4855 

83.2 

-15.37 

34.46 

6.49 

11.28 

3.49 

34.32 

> 

8.14 

< 

7.09 

< 

-2.09 

Q4+0.19 

5.43 (6) 

J0205+6449 

65.7 

-12.72 

37.41 

3.74 

12.56 

5.06 

34.39 

> 

8.50 

< 

6.58 

< 

— 2.60 

2Q q. 5+0.67 
2y.ycJ_o 33 

11.47 (7) 

J0248+6021 

217.1 

-13.26 

35.32 

4.80 

12.54 

3.50 

34.39 

> 

7.69 

< 

6.99 

< 

— 2.19 

31 filt'O-OS 

1.60 (6) 

J0357+3205 

444.1 

-13.88 

33.77 

5.73 

12.38 

2.41 

<35.71 

6.92 

+ 0.01 
-0.14 

7.93 

+ 0.30 
-0.02 

-1. 

+ 0.30 
^'^-0.02 

33 14+0 01 

oo. j.-»_0 19 

3.96 (7) 

J0534+2200 

33.6 

-12.38 

38.64 

3.10 

12.58 

5.96 

35.79 

> 

8.71 

< 

6.69 

< 

— 2.49 

31 fi«^+0.01 

•^l-bb-o.og 

32.53 (9) 

30631+1036 

287.8 

-12.98 

35.23 

4.64 

12.74 

3.33 

33.75 

> 

8.47 

< 

6.19 

< 

-2.99 


5.35 (7) 

J0633+0632 

297.4 

-13.10 

35.08 

4,77 

12.69 

3.24 

<35.93 

7.72 

+ 0.22 
-0.14 

7.05 

+ 0.22 
-0.30 

-2. 

10 + 0.22 
^^-0.30 

32.83ti;:;i 

7.66 (7) 

J0633+1746 

237.1 

-13.96 

34.20 

5.53 

12.20 

3.05 

34.50 

7.66 

-0.07 

7.15 

-0.01 

-1. 

Q0 + ^' -*-b 
^'^-O.Ol 

31 

Q4 

194.30 (11) 

J0659+1414 

384.9 

-13.26 

34.58 

5.04 

12.67 

2.88 

32.37 

6.42 

-0.42 

8.06 

-0.40 

-1. 

"^^-0.40 

30.9+“:™ 

3.68 (5) 

J0734-1559 

155.1 

-13.90 

35.11 

5.29 

12.15 

3.54 

<35.85 

> 

7.60 

< 

7.18 

< 

-2.00 

32.95j;0;;2 

19.25 (8) 

J0835-4510 

89.4 

-12.90 

36.84 

4.05 

12.53 

4.64 

34.95 

> 

8.03 

< 

7.19 

< 

-1.99 

31.19+0:5; 

1418.00 (11) 

J0908-4913 

106.8 

-13.82 

35.69 

5.05 

12.11 

3.99 

34.54 

7.22 

+ 0.23 
-0.19 

7.31 

+ 0.33 
-0.30 

-1. 

40 + 0.33 
^^-0.30 

3i.8oj:0;“ 

1.72 (6) 

J1016-5857 

107.4 

-13.09 

36.41 

4.32 

12.48 

4.34 

34.74 

> 

7.90 

< 

7.40 

< 

— 1.78 

30.86^0;^ 

3.68 (6) 

J1023-5746 

111.5 

-12.42 

37.04 

3.66 

12.82 

4.64 

<36.82 

> 

8.94 

< 

6.08 

< 

-3.10 

39 oq+0.30 
"‘^•'^^-0.23 

6.06 (9) 

J1028-5819 

91.4 

-13.79 

35.92 

4.95 

12.08 

4.17 

35.20 

> 

8.45 

< 

6.69 

< 

— 2.49 

31 30'^^'^^ 

oj..ou_o 08 

9.05 (10) 

J1044-5737 

139.0 

-13.25 

35.92 

4.59 

12.45 

3.99 

<36.74 

> 

7.93 

< 

6.99 

< 

— 2.19 

32 Qfi+0.48 

d2.yb_o 03 

18.97 (8) 

J1048-5832 

123.7 

-13.02 

36.30 

4.31 

12.54 

4.23 

35.25 

> 

7.91 

< 

7.19 

< 

-1.99 

ji.t)U_o 01 

12.09 (10) 

J1057-5226 

197.1 

-14.24 

34.48 

5.73 

12.04 

3.11 

33.64 

7.40 

oo 
66 
+ 1 

7.47 

oo 
66 
+ 1 

-1. 

o 

1 + 
o o 
66 

OU.^J_0 01 

5.89 (8) 

J1105-6107 

63.2 

-13.80 

36.40 

4.80 

12.00 

4.57 

35.16 

> 

8.26 

< 

6.90 

< 

— 2.28 

3 1 bH'TU • 
01-28-0.54 

2.08 (6) 

J1119-6127 

408.7 

-11.39 

36.36 

3.21 

13.61 

3.74 

35.78 

> 

8.31 

< 

6.29 

< 

-2.89 

32 04"^^'^^ 
o^.u4_o oe 

5.31 (8) 

J1124-5916 

135.5 

-12.12 

37.08 

3.46 

13.00 

4.58 

35.23 

> 

9.32 

< 

5.58 

< 

— 3.60 

30 48+0.iiU 

OU.4B_o 02 

4.73 (8) 

J1135-6055 

114.5 

-13.11 

36.32 

4.36 

12.48 

4.27 

<36.29 

> 

8.32 

< 

6.68 

< 

-2.50 

32.38j;0;2^ 

1.98 (6) 

J1413-6205 

109.7 

-13.56 

35.91 

4.80 

12.26 

4.09 

<36.94 

> 

8.06 

< 

7.20 

< 

-1.98 

33.0810;“; 

5.09 (8) 

J1418-6058 

110.6 

-12.77 

36.69 

4.02 

12.64 

4.47 

34.97 

> 

8.12 

< 

6.99 

< 

— 2.19 

3 1 93 + 0 • U / 

4.88 (8) 

J1420-6048 

68.2 

-13.08 

37.00 

4.11 

12.38 

4.84 

35.81 

> 

8.76 

< 

6.39 

< 

— 2.79 

31 9S + ^'^1 
^^•^®-0.39 

4.04 (8) 

J1429-5911 

115.8 

-13.52 

35.89 

4.78 

12.28 

4.05 

<36.66 

> 

7.76 

< 

7.29 

< 

-1.89 

33 98+ 0-07 

oo.^o_0_08 

6.55 (7) 

J1459-6053 

103.2 

-13.60 

35.96 

4.81 

12.20 

4.14 

<36.88 

> 

8.76 

< 

6.18 

< 

-3.00 

33 00~^0-i4 
jj.uu_o 03 

4.89 (8) 

J1509-5850 

88.9 

-14.04 

35.72 

5.18 

11.96 

4.08 

35.02 

> 

8.05 

< 

7.01 

< 

— 2.17 

o+oolSII 

45.43 (8) 

J1620-4927 

171,9 

-13.98 

34.91 

5.41 

12.15 

3.39 

<37.04 

> 

7.82 

< 

7.19 

< 

-1.99 

33.83l0:;0 

9.44 (9) 

J1709-4429 

102.5 

-13.03 

36.53 

4.24 

12.49 

4.43 

35.93 

> 

8.10 

< 

7.09 

< 

-2.09 

32.0912:55 

135.00 (11) 

J1718-3825 

74.7 

-13.88 

36.08 

4.95 

12.00 

4.35 

35.14 

> 

7.65 

< 

7.49 

< 

— 1.69 

O7+0.10 

di.87_o 70 

4,62 (6) 

J1732-3131 

196.5 

-13.55 

35.18 

5.05 

12.38 

3.46 

33.94 

7.56 

-0.04 

7.47 

-|-U • J- i 
-0.01 

-1. 

^^-0.01 

3(1 oa + O.Ul 

30.28-005 

0.80 (7) 

J1741-2054 

413.7 

-13.77 

33.98 

5.59 

12.43 

2.54 

33.31 

6.92 

+ 0.10 
-0.10 

7.99 

+ 0.23 
-0.20 

-1. 

+ 0.23 
■^'^-0.20 

54+0.15 
OO.t>4_0 13 

7.42 (7) 

J1746-3239 

199.5 

-14.18 

34.52 

5.68 

12.08 

3.13 

<36.74 

> 

7.53 

< 

7.29 

< 

-1.89 

34.08l0:g 

6.38 (6) 

J1747-2958 

98.8 

-13.21 

36.40 

4.41 

12.40 

4.38 

35.76 

> 

8.02 

< 

6.98 

< 

— 2.20 

31.7810;“ 

4.73 (8) 

J1803-2149 

106.3 

-13.71 

35.81 

4.94 

12.18 

4.05 

<36.85 

> 

8.27 

< 

6.79 

< 

-2.39 

33.09i0:;5 

6.49 (8) 

J1809-2332 

146.8 

-13.46 

35.63 

4.83 

12.36 

3.82 

35.21 

> 

7.94 

< 

7.08 

< 

— 2.10 

31.74l0:0; 

11.58 (9) 

J1813-1246 

48.1 

-13.75 

36.79 

4.64 

11.97 

4.89 

<37.27 

> 

9.22 

< 

5.98 

< 

-3.20 

32.6812:52 

8.11 (8) 

J1826-1256 

110.2 

-12.92 

36.56 

4.16 

12.57 

4.40 

<37.45 

> 

7.73 

< 

7.39 

< 

— 1.79 

oj.yB_o_o8 

17.43 (8) 

J1833-1034 

61.9 

-12.69 

37.53 

3.69 

12.56 

5.14 

35.19 

> 

7.80 

< 

7.39 

< 

— 1.79 

31.9212:52 

1.57 (6) 

J1836+5925 

173.3 

-14.82 

34.04 

6.26 

11.72 

2.96 

34.31 

> 

7.70 

< 

7.29 

< 

-1.89 

31.0912:2; 

22.09 (9) 

J1838-0537 

145.7 

-12.33 

36.77 

3.70 

12.92 

4.39 

<37.11 

> 

8.77 

< 

6.28 

< 

— 2.90 

55+0.40 

04i.t>t>_o_o6 

8.43 (8) 

J1907+0602 

106.6 

-13.06 

36.45 

4.29 

12.49 

4.37 

35.50 

> 

8.00 

< 

7.09 

< 

— 2.09 

31 1 R + 0 • ^ ^ 

31-18-0.03 

24.13 (9) 

J1952+3252 

39.5 

-14.24 

36.57 

5.03 

11.68 

4.86 

34.82 

> 

8.07 

< 

7.39 

< 

— 1.79 

30.9312:22 

12.61 (9) 

J1954+2836 

92.7 

-13.67 

36.04 

4.84 

12.15 

4.21 

<36.63 

> 

8.10 

< 

7.08 

< 

— 2.10 

32.8512:25 

6.57 (9) 

J1957+5033 

374.8 

-14.17 

33.71 

5.94 

12.20 

2.45 

<35.82 

7.06 

1 + 
o o 

7.65 

+ 0.33 
-0.31 

-1. 

+ 0.33 
'^'^-0.31 

33 3Q+0-21 

BB.By_0 28 

6.25 (7) 

J1958+2846 

290.4 

-12.67 

35.53 

4.34 

12.90 

3.47 

<36.57 

> 

7.29 

< 

7.70 

< 

-1.48 

33 47”^0 - 0 (’ 
_o.i9 

20.68 (8) 

J2021 + 3651 

103.7 

-13.02 

36.53 

4,23 

12.51 

4.42 

36.77 

> 

7.94 

< 

7.19 

< 

-1.99 

33 1 5+0.08 

35-lb-o.oi 

24.60 (10) 

J2021 + 4026 

265.3 

-13.27 

35.04 

4.89 

12.58 

3.28 

35.41 

> 

7.65 

< 

7.09 

< 

— 2.09 


58.07 (11) 

J2028+3332 

176.7 

-14.31 

34.54 

5.76 

11.97 

3.20 

<36.32 

7.70 

+ 0.21 
-0.21 

7.03 

+ 0.43 
-0.30 

-1. 

90 + 0.43 
^ —0.30 

33.1312:;; 

8.56 (8) 

J2030+3641 

200.1 

-14.19 

34.51 

5.69 

12.08 

3.12 

34.53 

7.46 

+ 0.26 
-0.27 

7.68 

+ 2.30 
-0.61 

-1. 

on + 2.30 

30.8812:25 

5.61 (7) 

J2030+4415 

227.1 

-14.19 

34.34 

5.74 

12.08 

2.98 

<36.23 

> 

7.73 

< 

6.88 

< 

-2.30 

33 QR ^ ^ 

BJ.yB_o 10 

1.33 (6) 

J2032+4127 

143.2 

-13.69 

35.43 

5.05 

12.23 

3.73 

35.23 

8.00 

+ 0.12 
-0.18 

6.93 

+ 0.41 
-0.20 

-1. 

90 + 0.41 
^'^-0.20 

31 34+0-12 

b^-'54_o 22 

6.15 (8) 

J2055+2539 

319.6 

-14.39 

33.70 

6.09 

12.08 

3.51 

<36.18 

7.10 

-0.05 

7.88 

-0.21 

-1. 

30 + 0'11 
’^^-0.21 

OOl+S'.TO 

3.43 (7) 

J2111 + 4606 

157.8 

-12.84 

36.15 

4.24 

12.68 

4.05 

<36.06 

> 

8.05 

< 

6.99 

< 

— 2.19 

32.26i2:;2 

9.10 (9) 

J2139+4716 

282.8 

-14.74 

33.49 

6.40 

11.86 

2.47 

<35.75 

> 

7.15 

< 

7.71 

< 

— 1.47 


6.83 (6) 

J2229+6114 

51.6 

-13.11 

37.34 

4.02 

12.30 

5.13 

34,29 

9.22 

+ 0.19 
-0.60 

5.29 

+ 0.60 
-0.20 

-3. 

10 + 0.60 
^^-0.20 

29.6512:52 

13.01 (10) 

J2238+3903 

162.7 

-13.01 

35.95 

4.42 

12.60 

3.93 

<36.07 

> 

7.67 

< 

7.29 

< 

-1.89 

32.7912:52 

7,75 (9) 
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Table 3. Same as Table [5] but for the 22 MSPs of our sample. 


Pulsar 

P 

log P 

log E 

log T 

log Bs 

log Bi^ 

log L 

log P|| 

log a;o 


log A^o 

^min^-^bins) 


|ms| 


[erg/s] 

[y] 

[G] 

[G] 

|.rg/B] 

[V/m] 

[cm] 



J0030 + 0451 

4.9 

-19.99 

33.53 

9.88 

8.36 

4.25 

32.76 

0 eR”l"0.14 

= =7+0.22 
^•°^-0.21 

-i.80l«;== 

28.53l°;j3 

10.78 (8) 

J0034-0534 

1.9 

-20.70 

34.46 

10.18 

8.00 

5.12 

32.75 

> 9.01 

< 5.38 

< -1.99 

28.38+“'“ 

1.45 (6) 

JOlOl-6422 

2.6 

-20.32 

34.04 

9.93 

8.04 

4.77 

32.58 

Q yp-f-0.19 
®-^°-0.20 

. qq+2.10 

^•yy_o.5i 

10+2-10 

27 51+0-81 
27.51_o 16 

1.99 (6) 

J0218 + 4232 

2.3 

-19.11 

35.40 

8.67 

8.63 

5.51 

34.58 

> 9.84 

< 4.47 

< -2.90 

29.5i+°n5 

1.37 (9) 

J0340+4130 

3.3 

-20.23 

33.81 

9.95 

8.15 

4.56 

33.86 

„ R4+0.42 

®-®^-0.37 

60+2-60 

^•^^-0.70 

-1 60+2-00 
-^•'^'-^-0.70 

97+0.61 

6.12 (6) 

J0437-4715 

5.8 

-19.27 

34.04 

9.23 

8.75 

4.42 

31.69 

> 8.19 

< 5.88 

< -1.49 

28.06+“'j® 

4.04 (7) 

J0613-0200 

3.1 

-20.02 

34.11 

9.71 

8.26 

4.73 

33.46 


5 27+°-^^ 
^•2^-0.50 

-1 90+0-44 

-^-y^-0.50 

OQ Q4+0.26 

28.84_q 28 

4.99 (6) 

J0614-3329 

3.2 

-19.75 

34.32 

9.45 

8.38 

4.83 

34.67 

> 9.04 

<5.47 

< -1.90 

29'86to,S| 

12.62 (10) 

J0751 + 1807 

3.5 

-20.11 

33.86 

9.85 

8.23 

4.56 

32.40 

> 9.00 

< 5.39 

< -1.98 

28.02+0;9^ 

2.22 (6) 

J1124-3653 

2.4 

-20.24 

34.23 

9.82 

8.08 

4.90 

33.63 

> 9.02 

< 5.59 

< -1.78 

98-95 + 0 27 

1.72 (6) 

J1231-1411 

3.7 

-19.67 

34.20 

9.44 

8.45 

4.71 

33.37 

Q qo+O-OS 
®-^^-0.12 

. .,. + 0.20 
'^•'^°-0.11 

— 1 90+0-20 

28.76l“;5| 

12.81 (9) 

J1514-4946 

3.6 

-19.73 

34.20 

9.48 

8.41 

4.71 

33.68 

> 9.20 

< 5.18 

< -2.19 

29.05+9-™ 

1.85 (8) 

J1614.2230 

3.2 

-20.02 

34.08 

9.72 

8.26 

4.70 

33.09 

Q CQ+O-IS 

»-&®_0.14 

6 oq + 1.80 

6-6»_o.73 

n sn+l-OO 

-0.80_o 73 

2T.9e±°ll 

3.29 (6) 

J1744-1134 

4.1 

-20.05 

33.71 

9.86 

8.28 

4.41 

32.83 

> 8.42 

< 5.78 

< -1.59 

29.00 

2.06 (6) 

J1810+1744 

1.7 

-20.34 

34.57 

9.77 

7.95 

5.23 

34.05 

> 9.41 

< 4.87 

< -2.50 

29 68+'-^'^ ' 
^y.D5_o 04 

2.28 (7) 

J1902-5105 

1.7 

-20.05 

34.86 

9.48 

8.08 

5.37 

33.56 

> 9.52 

< 4.87 

< -2.50 

OQ Bfi+0-04 

2»-»6 _o.o 8 

2.54 (6) 

J2017+0603 

2.9 

-20.08 

34.11 

9.74 

8.20 

4.77 

33.99 

9 22+°-^® 

y-^^_o.i5 

5 14+0-33 
^■•'•^-0.30 

_ 00+0-03 

.^.uu_o 30 

28 96+0-1"^ 
^».yD_0 20 

4.26 (8) 

J2043+1711 

2.4 

-20.24 

34.20 

9.82 

8.08 

4.90 

34.00 

> 9.13 

< 5.28 

< -2.09 

29.47+0-14 

4.60 (8) 

J2124-3358 

4.9 

-19.69 

33.84 

9.58 

8.51 

4.40 

32.60 

Q cn + 0.08 

6 37+0.61 
‘"•‘^^-0.31 

— 1 00+0-01 
-0.31 

27 65+0-12 
27-b5_o 10 

2.07 (7) 

J2214+3000 

3.1 

-19.82 

34.30 

9.51 

8.34 

4.83 

33.97 

0 qq-|-0.09 
^■^^- 0.11 

= 37+0.21 
'^•■^^-0.11 

sn+0-21 

-1.80_o.ii 

29 33+0-06 
zy.oo_o 11 

3.72 (7) 

J2241-5236 

2.2 

-20.06 

34.51 

9.60 

8.15 

5.08 

33.02 

> 9.04 

< 5.58 

< -1.79 

28.15t“;“ 

2.09 (7) 

J2302 + 4442 

5.2 

-19.88 

33.57 

9.79 

8.43 

4.24 

33.79 

„ -n+0.11 

®-”°-0.14 

5 79+0-51 
*^^-0.21 

_ 60+0-0^ 
i.ou_o 21 

29.08 

4.48 (8) 


the boundaries Xin and Xout, and the parameters describ¬ 
ing the decreasing of B, and the increasing of along the 
trajectory of the particles. The relative irrelevance of these 
parameters is, in itself, an interesting result. On the one 
hand, it makes the model more robust. On the other hand, 
it tells that the spectra are not particularly sensitive to de¬ 
tails like the precise form of the magnetic field decay in the 
magnetosphere, that would seem otherwise key. 


3 FITTING THE PULSAR SPECTRA 
3.1 Sample definition and fitting 

We consider the publicly available, Fermi-LAT processed 
data for the phase-averaged pulsar spectra contained in the 
2PC. As a good-quality criterion, and not to consider less 
constraining data, we select those sources having measured 
flux (not upper limits) in at least five consecutive energy 
bins. Such selection criterion somewhat favors the YPs, since 
they are brighter in 7 -rays. Among the 81 sources of our 
sample (out of the 117 of the 2PC), 59 are YPs (Table [21) 
and 22 are MSPs (Table (3]). 

We apply our models exactly in the same way to each 
YP and MSP in the sample. We neither make any distinc¬ 
tion between radio-quiet or radio-loud. Thus, we constrain 
the values of the three relevant parameters, Fy, xo, and No, 
by comparing the SC radiation expected from our models 
once they assume a particular value with the observational 
data. The fitting procedure is described in detail in Ap- 
pendixj^ We span a range of several orders of magnitude in 
the values of the three mentioned parameters, thus consid¬ 
ering a grid of thousands of models for every pulsar before 
selecting the best fitting set of parameters. For each model 
tested, we compute the particle trajectory in the gap and 
obtain the gap-integrated spectra according to Eq. (|T]). To 


quantitatively compare model and data we compute (see 
Appendix m for details). 

We model only the particles directed to the observer, 
in an effective way, as stated in §2 with No being the nor¬ 
malization of the number of particles that are moving in the 
direction of the observer. The fits show what collection of 
the three model parameters can fit the data. The gap width 
is not a parameter in our model, since we have no geometry 
included. 

3.2 Fitting results and model assessment 

Tables [2] and [3] show the timing and spectral properties of 
the Fermi-LAT pulsars in our sample, together with the 
best-fitting parameters inferred by applying our model. Ap- 
pendixlBlpresents the case by case results. Figs. IB lllBSI show 
the observed E^dN/dE spectra (where dN/dE is the spec¬ 
tral photon distribution radiated by the source) for each of 
the 81 pulsars in the sample together with the corresponding 
best-fitting models. The theoretical spectra for YPs (MSPs) 
are plotted there in red (blue), a color code we maintain in 
all sections. 

The best-fitting models show quantitative agreement 
with data. With just three parameters, the model copes rel¬ 
atively well with data of dozens of pulsars, disregarding their 
type (YP or MSP) and spanning several orders of magnitude 
in any of their main characteristics. 

3.2.1 Analysis of the residuals 

The highest values of x^ occur for a few, very bright YPs: 
J0835-4510 (Vela pulsar, > 1000, dof = 8 ), J0633-fl746 
(Geminga = 194.3, dof = 8 ), and J1709-4429 (x^ = 135, 
dof = 8 ). They show residuals at high energies, similarly to 
the Crab pulsar (x^ = 32.53, dof — 6 ). There are two main 
reasons producing the large values of x^ in these cases. On 








Fits of pulsar 'j-ray spectra 5 



Travelled distonce/Ri^ 



11 


9 



OJ 

o 

o 

Q. 



10 “^ 10 “^ 10 “^ 10 “^ 10 “’ 10 ° 
Travelled distance/Ri^ 


10”'* lO"-" 10"^ 10” 

Travelled distanae/R|g 


11 




11 


9 



OJ 

o 

o 

Q. 


Figure 1. Evolution of the trajectory parameters for YPs (left-hand panels) and MSPs (right-hand): Lorentz factor of the particles in the 
gap, and g value, determining the dominance of synchrotron-like (g 1 ) or curvature-like (g <g 1 ) radiation along the traveled distance 
in the gap (x), normalized with the corresponding light cylinder radius. The color scheme corresponds to the value of the parallel electric 
field. 


the one hand, being these the brightest pulsars, the rela¬ 
tive error bars at the low-energy part of the spectra are 
extremely small (better than 1 %), and this entails a rela¬ 
tively large deviation from theoretical curves even when the 
overall description as shown in the corresponding spectral 
plot of Appendix m is appropriate. On the other hand, the 
high-energy end of the spectra of these pulsars is decaying 
slower than what the model predicts (Fermi-LAT data is 
better described with a power-law with a sub-exponential 
cutoff, dubbed PLEC in the 2PC). This may, in fact, be a 
spurious effect of considering phase-averaged spectra. 

The latter is the case for Geminga, where the set of 
fine-binned phase-resolved spectra has shown the absenc e 
of sub-exponential decays at all phases llAbdo et al.ll 2010 ^ ). 
i.e., being described by power-law with exponential cut-off 
models (PLECl in the 2PC). The phase-averaged, harder 
Geminga spectrum arises as a consequence of summing up 
phase-resolved contributi ons that can be indiv idually well 
described by a SG model dVigano fc Torre j|2015l ). Two more 
cases are similar to Geminga, with which they share simi¬ 
lar values of P and P: J0007+7303 in CTAl = 41.37, 
dof = 8 ) and J2021+4026 (y^ = 58.07, dof = 8 ). Eor both, 
the pha se-resolved spectra can always be fitte d by PLECl 
models llAbdo et al.lf20l3 : lAllafort et al.ll2013l ). 

Crab and Vela are different, how ever, since they 
show pulsed emis sion above tens of GeV (lAliu et al.ll2008l : 
iLeung et al.ll2014l ). and, in at least one phase, their phase- 
resolved spectra is harder than an exponentially decaying 
one. Such high-energy tails are very likely a signature of in¬ 


verse Compton processes contributing to the total yield and 
cannot be described well by just a SC model, particularly 
at high energies. We have discussed the cases of Geminga, 
Crab, and Vela pulsars in more d etail using phase-resolved 
spectra in IVigand fc TorresI (l2015l) . 

In some cases, the first couple of data points show 
non-negligible residuals within a few ct, or upper limits 
apparently incompatible with our model. Eor instance, in 
J1907+0602 and J1958+2846, the low-energy slope is appar¬ 
ently harder than what the SC radiation of the best-fitting 
model produces. If these discontinuous points were real (e.g., 
the first point in J1509-5850), then the phenomenology of 7 - 
ray emission is more complicated than what this SC model is 
able to describe, and perhaps include more than one compo¬ 
nent. We note, however, that diffuse model uncertainties can 
especially affect the first couple of low-energy data points. 
In this regard, we hope that pulsar results using Pass 8 in 
analyzing Fermi-LAT data will bring additional insight. 


3.2.2 Trajectories and spectra 

It is interesting to relate the best-fitting spectra with the 
corresponding particle dynamics along the gap, shown in 
Fig. [T] for the YPs (left) and the MSPs (right). We show 
the evolution of their Lorentz factor, P (top), and their SC 
parameter, ^ (bottom) -see Eq. IA 8 I and Appendix for 
details-, as a function of travelled distance along the field 
line. Colors indicate the intensity of the accelerating field: 
the larger Fy, the larger the final T. 
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We see the overall similarities of models, despite the dif¬ 
ferent spectra they produce. Particle trajectories start being 
dominated by the loss of perpendicular momentum (initial 
plateaux in P, ^ 3 > 1 ), then the electric acceleration boost 
them up, till they reach a balance between the electric accel¬ 
eration, and the SC losses, dominated by the loss of longitu¬ 
dinal momentum (curvature regime, ^ <C 1). The transition 
between the different ^-regimes happens at about ( 10 ~®- 
10~^)i?ic. The larger is the value of the accelerating electric 
field, the sooner this transition happens. At this same (cor¬ 
responding) location (i.e., at physically different location for 
each of the pulsars, since it depends on the corresponding 
Ric) the Lorentz factor saturates. 

In MSPs, despite having a larger best-fitting value of 
, particles reach lower (or similar) values of P, because the 
radiative losses are stronger, being the light cy linder smaller 
(and, consequently also their curvature radius. Ivigand et al.l 
l2015bl l. The different saturated values of the Lorentz factors 
between MSPs and YPs forms a continuing set similar to iJ|| 
values, discussed in the next section. 

If the effective particle distribution were uniform, parti- 
cles with saturated val ues of T would dominate the spectra 
dVigano fc Torr^l2015l L since their radiation is more intense 
(higher flux) and harder (larger energy peak). However, this 
could not explain the variety of spectral slopes. The flat¬ 
ness of spectra require that less energetic particles, having 
lower P and larger to significantly contribute as well. In 
our model, the parameter xq regulates the relative weight of 
this contribution. 

3.2.3 Error estimate and degeneracy 

Figs. IBhllBlOl of Appendix iBlshow the contours of X^/Xmin 'n 
the plane log i5||-log(a;o/i?ic). Again, red and blue labels help 
distinguish the pulsar class. The contour plots measure the 
deviation from the minimum attainable value of in the 
model grid, i.e., when the parameters differ from the best¬ 
fitting set. In each of the panels, we mark with a white cross 
the location of the best-fitting model (i.e., where = Xmin)- 

In order to estimate the errors on parameters, we con¬ 
sider the region of the contour plots where x^/Xmin ^ 1-3 
(approximatively, at the boundary between black and red re¬ 
gions). Then, we take the minima and maxima of the three 
parameters within the enclosed region to estimate the er¬ 
rors. For pulsars having a well-defined region, typical rela¬ 
tive errors on the inferred best-value parameters of about 
5Zo(/(iJ||)/log(i7||) ~ 5 — 10%, ^log(xo)/log(xo) ~ 30 — 50%, 
and (5 log(Ao)/log(Ao) ~ 1%- These ranges of uncertainties 
are acceptable since we explore several orders of magnitude 
in the three parameters for each of the pulsars. 14 YPs and 
10 MSPs have well-constrained parameters, and thus it is 
possible to study trends in the population using just these 
(filled circles In all the plots). 

Many pulsars show a very elongated (or linear) shape 
in their corresponding (log-log) contour plot of Figs. IB 6 I 
IBIOI pointing to a clear anti-correlation between Sy and xq. 
In these cases, there is a systematic degeneracy in choosing 
the best-fitting model. This can be understood when one 
considers that at a larger E\\, particles are pushed to larger 
energies; while on the contrary, a smaller xo tends to neglect 
the part of the trajectories where particles attain the largest 
energies weighting more heavily on the synchrotron domi¬ 


nated part and thus compensating for the previous effect. 
This anti-correlation, Sy ~ for large enough E^\, can 
be analytically proven to hold from the equation of motions 
themselves (see Appendix|^for an approximate derivation), 
when the electric term dominates over the SC loss term. It 
does not hold whenever particles with larger F are dominat¬ 
ing the spectrum (larger xo): in these cases, the degeneracy 
is partly broken and the value of By is more constrained 
(see the contours with more vertically-oriented shapes, e.g. 
JOlOl-6422). 

From the trajectories shown in Fig. [T] one can under¬ 
stand that a larger By displaces the r(x) curves to the left; 
thus, to reach the same observed spectrum, one has to com¬ 
pensate by convolving a particle distribution weighted to¬ 
wards lower values of x. In such cases, the exact location of 
the minimum found by our pipeline is slightly dependent on 
the resolution and range of the explored grid of parameters: 
all the pairs of parameters inside the black and red valley In 
the pulsars logi?y-log(xo/i?ic) contour plots of Appendix [Bl 
are closely as good as the ones marked by the white cross. 

The degeneracy can be partially solved on physical 
grounds. In particular, we think that within the constrained 
range, the lowest values of the allowed range of i?y and the 
corresponding largest values of xo are more physically mean¬ 
ingful. As a matter of fact, it is difficult to justify large val¬ 
ues of E\\ exceeding the local value of magnetic field, B\c 
by orders of magnitude, and, at the same times, values of 
Xo of the order of meters (with an accordingly very large 
particle density). In Tables [2] and [3l the pulsars presenting 
a stronger degree of degeneracy in the selection of the best¬ 
fitting values (over two thirds of the sample) are represented 
by lower (upper) limit on By (xo). Such lower/upper limits 
are also employed in the correlation plots hereafter, identi¬ 
fied by empty circles. 


4 CORRELATIONS AND TRENDS 

In this section we explore correlations between the best¬ 
fitting model parameters and the timing and timing-derived 
properties of pulsars, the best phenomenological fit (PLECl 
model values) as quoted in the 2PC, and (briefly) the pulse 
profile characteristics. 

4.1 Values of parameters and self-correlations 

Fig. [3 summarizes the distributions of our best-fitting SC 
model parameters. As Fig. [5] and the Tables [5] and [3] show, 
MSPs can typically have up to one order of magnitude larger 
Ey when compared with YPs. Large electric fields, Ey ~ 
10 ®-® — 10^° V/m (compared with Ey ~ 10®'® — 10® V/m 
for YPs) are required in MSPs in order to compensate for 
their much smaller radius of curvature (due to the much 
closer light cylinder) and let particles be energetic enough 
(r ~ 10^-10®, see Fig. [1} as to be able to produce ~ GeV 
photons. The lowest (highest) values of Ey are found for YPs 
(MSPs), respectively. 

Typical values of xo /Ric are consistently low across both 
sub-samples of pulsars, (or, equivalently, xq ~ 1 - 

1000 km), as wa s already the case fo r Crab, Geminga, 
and Vela pulsars llVigano fc Torre^l2015l L These values are 
needed to accommodate the many, soft spectral indices of 
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Figure 2. Distributions and correlations between our best-fitting SC model parameters. Hereafter, YPs and MSPs are marked by red 
and blue points, respectively, and empty circles represent the cases where i?|| and xq can only be constrained by lower and upper limits, 
respectively. Therefore, empty points could actually be moved along i?|| ~ (l/xg) towards larger (smaller) values of i?|| (xq). See text 
for a detailed discussion. Lines represent linear fits (shown only when the corresponding Pearson coefficient is r > 0.85, in black for the 
whole sample, in red (blue) for YPs (MSPs)) to the non-degenerate pulsar parameters only (filled circles). 


the 7 -ray pulsars below the cutoff energies, as observed by 
Fermi-LAT. They indicate a common feature: the 7 -ray ra¬ 
diation is not produced only by particles with saturated val¬ 
ues of Lorentz factor. Instead, most photons need to be pro¬ 
duced while particles have non-negligible perpendicular mo¬ 
mentum, thus radiating mostly synchrotron emission. The 
lowest (highest) values of xq are found for MSPs (YPs), re¬ 
spectively. 

The values of No clearly separate YPs from MSPs. This 
is a reflection of the MSPs being less luminous in general, 
which directly implies that all MSPs have smaller A^o (in the 
extreme case, for up to six orders of magnitude). Physically, 
this could be related either to the much reduced magne- 
tospheric volume in MSPs, and/or to the lower rotational 
energy losses, which represent the ultimate energetic budget 
partially converted in particle acceleration, and, therefore, 
in visible radiation. 

We note that, if the number of effective particles we in¬ 
fer were all born at the inner boundary of the gap, then the 
inferred flux of particle, (A^o/ 2 ;o)c, is much larger (between 
a factor of a few up to 4 orders of magnitude) than the esti¬ 
mated Goldreich-Julian flux coming from t he open field line 
regio n and hitting the polar cap (see e.g. jZhang fc Chena 
Il997t) ). Our larger number of particles could be an indica¬ 
tion of the oversimplification of our model (e.g., particles are 
generated in different altitudes), or that the volume occu¬ 
pied by the gap is much larger than thought, or simply that 
cascading is very active (large multiplicity factor, similarly 


to pulsar wind neb ulae models, e.g.. lBucciantini et al.|[201ll : 
[Torres et al.ll2014l) . Only numerical simulations considering 
all the magnetosphere could provide detailed answers. How 
active can the cascading be before starting to significantly 
screen the accelerating electric field also needs simulations 
to tell. 

Among the three parameters, we find that only i?|| and 
xo maintain a significant anti-correlation, which is sepa¬ 
rately visible and continuing across both sub-samples (upper 
left-hand panel of Fig. n. This is confirmed by the Pearson 
coefficient shown in Tabled] Note that the correlation lines 
in Fig. [ 2 ] and all subsequent ones, as well as the correlation 
parameters in Table [4] are obtained considering the filled 
circles only, i.e., the pulsars for which all F/n, xo, and No are 
well-constrained. The degenerate pulsar cases match these 
correlations very well in all cases. 

The global correlation in the upper left-hand panel of 
Fig. [2]shows all the best-fitting solutions together. It should 
not be confused with the individual correlations found in 
the same plane for many pulsars in Figs. lB6IB10l Those con¬ 
tours cover different ranges of parameters from pulsar to pul¬ 
sar (which is natural since each pulsar has different timing 
properties), but the linear parts of the contours have similar 
slopes, as already discussed in il3.2.3l and in Appendix 
Note that the global trend will be maintained despite the 
degeneracy in the values of parameters in many individual 
cases, since the latter would move the points along a similar 
diagonal line in the global plane log F/y-log xo of Fig. [2] In 
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Figure 3. Distribution of the three parameters fitted in our models as a function of P (top panels) and P (bottom). Color and empty/filled 
codings are the same as Fig. [2] 


other words, the particular degeneracy in determining the 
specific pulsar best-fitting model does not affect the overall 
trend. Since we find no reason to expect any a-priori bias, 
this global trend is a novel result, which can highlight phys¬ 
ical implications. 

The correlation also points to the fact that particles 
are accelerated to similar Lorentz factors, as the rightmost 
part of Fig. 1 (top-left panel) show, between ~ 2 x 10^ and 
3 X 10®. These high values can be interpreted as those needed 
to start pair cascades, which ultimately support the gap. In 
fact, the product xqE\\ stays within 1 order of magnitude 
around 10^® V, a value that has been seen as the voltage 
that is limited by curvatur e radiation pair cascades, e.g., 
iHarding fc Musliinovl ll2002h . 

The linear anti-correlation is found also when xo is 
rescaled by i?ic, i?|| ~ {xo/Ric)~^ (see top right-hand panel), 
but, in this case, the two sub-samples get separated (due to 
the very different values of iiic). This reinforce the idea that 
the same underlying microphysical process is at work, being 
barely dependent on the macrophysical lengthscale of the 
fight cylinder. In other words, the observed spectra sistem- 
atically present features which can be explained by the local 
particle dynamics, more than by the global properties like 
the magnetospheric size (set by the spin period). 

4.2 Correlations with P and P 

Fig. [3] shows the distribution of the three parameters fitted 
in our models as a function of P and P. The two populations 
are clearly distinguishable. 

The distribution of the accelerating electric field is 
shown in the left panels. Taking the two sub-samples indi¬ 
vidually, i5|| hints to have increasing values for smaller spin 
periods. This is consistently reflected in the middle panel, 
where the hint is for a larger value of xo at larger periods. If 
true, such behavior can be understood as follows: the smaller 


is the period of the pulsar, the closer is the light cylinder to 
the surface of the neutron star (NS), and the smaller is the 
range in which the Lorentz factor needs to achieve values 
high enough as for particles to be able to emit 7 -rays with 
E ~ GeV. The two sub-samples are separated, indicating 
that the period alone is not the main parameter determin¬ 
ing the electric field (certain values of Sii can correspond to 
different values of P). No apparent correlation is seen be¬ 
tween i5|| or Xo and P. and no correlation between A^o and 
the measured timing properties is seen. A larger xo being 
obtained for a larger P probably reflects too the pair cas¬ 
cade physics: a smaller P have larger i?||. In our models, 
the latter requires a smaller xo to reach high voltages. Here 
there is the interesting possibility of an observational bias, 
since we have only Fermi-LAT survey data, the peak of the 
emission of all these pulsars are at ~l-few GeV, a larger 
Xo value for a similar i?|| would reduce the peak energy to 
mostly observationally-unexplored regimes in the MeV-few 
hundreds of keV region. 


4.3 Correlations with timing-derived properties 

Since YPs and MSPs span very different ranges of values 
of P and P, it is natural to search for possible correlations 
between the best-fitting parameters of the SC model and 
the timing properties of the pulsars. The panels of Fig. 0] 
explore these possible correlations against the characteristic 
age, the rotational energy loss, the surface magnetic field, 
and the magnetic field at the light cylinder, defined as 


Tc = 

r 

^ ’ 


(4) 

Erot 

= 3.9 X ergs"®. 

(5) 

Bs = 

= 6.4 X 10®® pi/2pi/2 

G, 

( 6 ) 

B\c = 

= 5.9 X io®p-5/2pi/2 

G, 

(7) 
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Figure 4. Search for correlations among the timing-derived parameters of the pulsars: r, Erot, Bs, and Bi,, with the three best-fitting SC 
model parameters. Color and empty/filled codings are the same as Fig. |2] The bold dashed lines represent the fits, when the correlation 
is significant. In the bottom left panel, the light dashed line identifies Ey = Bi^. The magnetic field in the accelerating region is even 
larger than this, and can be computed according to its defition in Table 1. See text for a discussion. 


where P is implicitly expressed in seconds and P in s s“^. 
Magnetic field values are estimated according to the widely- 
used vacuum spin-down formula, and assuming standard val¬ 
ues of moment of inertia I = 10^® g cm^, and NS radius 
i?* = 10 km. Since their actual physical values could vary, 
an intrinsic uncertainty by a factor of a few is expected on 
the inferred values Bs and Sic. Below we discuss these quan¬ 
tities and their correlations. All the estimates above are just 
different combination of powers of P and P. Other quanti¬ 
ties can be construed, like, e.g., the surface electric voltage 
AV^ = (Bs47r^iZ*)/(2cP^) but this, being proportional to 
p-3/2pi/2^ is directly proportional to E°ot- Thus, any cor¬ 


relation found with Erot would reflect in a similar one with 
the voltage. When we hnd a correlation of note, we list in 
Table 13] the parameters of a linear fitting, together with the 
Pearson test coefficient. 


4-3.1 Spin-down age Tc 

The range of values of Tc, Eq. m, spans more than seven 
orders of magnitude, well separating the two sub-samples: 
Tc ^ 10^-10® yr for the YPs, Tc ^ 10®-10^° yr for the 
MSPs. There is no visible correlation between Tc and any 
of the three model parameters, in any of the sub-samples. 
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We also recall that Tc is, more properly said, the spin-down 
timescale, and it can be a fair estimate of the real age only if 
the spin-down torque mechanism has not suffered variations 
with time (e.g., evolution of the large-scale dipolar mag¬ 
netic field, accretion phases...), and the present spin period 
is much larger than the period at birth. Note that the lat¬ 
ter assumption could not hold for some fast-spinning YPs 
and MSPs, which likely underwent accretion phases with 
different spin-down torque mechanisms. Thus, the lack of 
correlation is not a surprise. 


4-3.2 Spin-down energy loss E^ot 

The second row in Fig. [4] shows the search for correlations 
of the model parameters with the spin-down energy loss, 
Eq. ©. There appears a separation of the two sub-samples 
of pulsars, due to the discussed differences in iJ||, and to 
the fact that MSPs have lower spin-down power than most 
YPs. Values of Erot span more than five orders of magnitude, 
with Crab being by far the most energetic pulsar of the 
sample. Hints that for both populations, the larger is the 
spin-down, the larger is i?|| might be argued, but is not at 
the level of having the Pearson coefficient r > 0.85 for the 
well-constrained pulsars in both sub-samples. 


4-3.3 Surface magnetic field Bs 

The range of values of Bs, Eq. covers almost six orders 
of magnitude, again keeping the two sub-samples well sep¬ 
arated (third row in Fig. [J). MSPs are expected to be less 
magnetic, since they are much older than YPs, and their 
magnetic fields have been strongly dissipated. We find no 
correlation between Bs and any of the model parameters. 
This is expected if y-ray emission is coming from the outer 
magnetosphere. 


4-3-4 Magnetic field at light cylinder Hic 

The canonical value of B\c, Eq. 0 , defined simply by rescal¬ 
ing Bs, with (i?*/r)®, has to be taken as a rough estimate 
of the magnetic field in the gap. In fact, our model assumes 
a softer magne tic field decay (descr ibed by H = Bs{R-k/x)^, 
with b = 2.5, IVigano et alj l2015bl and Table [T]) , which is 
qualitatively more consistent with the numerical simulations 
of magnetospheres of rotating NSs. Nevertheless, the clas¬ 
sical estimate Sic still gives a relative idea of the magnetic 
field far from the surface across the sample, and we checked 
that all trends found with Bic would be found also assuming 
the scaling (i?*/r)^ ®. Moreover, in our models the spectra 
are quite insensitive to the exact value of b. 

Despite MSPs are less magnetized (smaller Bs), they 
have on average larger values of light cylinder field {Bw ~ 
10^-10® G) as compared to YPs (Bw ~ 10^-10® G), since 
their light cylinder is much closer, located at a few stellar 
radii. In the fourth row of Fig. [4] we explore correlations be¬ 
tween the model parameters and Bic- The uppermost outlier 
(in Sic) in these plots is the Grab pulsar. 

We find that the values of Bw correlate with the accel¬ 
erating field Fi|| resulting from our model fits. More interest¬ 
ingly, YPs and MSPs form an indistinct continuous in this 
plot (see Fig. 0. The inferred relation between Bic and Fi|| 


(quoted in Table 0 is very close to a direct proportionality, 
Fi|| oc Bic. This correlation, together with the fact that iJ|| 
does not correlate with Bs is an evidence that the accelerat¬ 
ing gap locates in the outer magnetosphere. However, note 
that since the correlation would appear also further away 
than the light cylinder, we cannot say whether the acceler¬ 
ating location is close to the light cylinder (as in the classical 
OG models), or farther away (as in the striped wind mod¬ 
els). 

Looking at correlations with the other two parameters, 
we consistently find a correlation such that an increasing 
value of xo implies a decreasing value of Bw- No correlation 
of Bic with No is visible. 

In Fig. 0 bottom left panel, we also show the line corre¬ 
sponding to having Sic = i?||. We have discussed our model 
compared wi t h thos e of the outer gap magnetosphere in 
IVigano et al.l (l2015bll (see espec i ally th e section 3.3 and Ta¬ 
ble 2 therein), and Vigand et al.l ll2015d L We discussed there 
in detail the many underlying assumptions of the outer gap, 
as well as the caveats in the computations of their main fea¬ 
tures. Among them, those applying to the formula quoted for 
iJ|l, which is used to define a death line separating gamma- 
ray emitting pulsars, i?|| oc Bwf^, where / is the fractional 
gap width, and / < 1. We concluded there that the specific 
values one can obtain for / from this formula has doubt¬ 
ful physical content. Sic is only an estimate of the value 
of the local magnetic field at light cylinder, assuming a de¬ 
pendence of a static dipole. However, neither the real radial 
dependence of the field nor the exact location of the acceler¬ 
ation place are known. In the outer gap models, in addition, 
there are a number of assumptions related with the photon 
fields that participate in the bootstrap to sustain the gap 
itself, and some problems related with the computation it¬ 
self (cross sections, resonant scattering, see the references 
quoted above for details), and all of that is carried over to 
the estimations of /. It is nevertheless interesting that the 
Hic vs. fit we provide in Table 0 and Fig. 0 is compara¬ 
ble to (although not strictly compatible with) a slope of 1 . 
However, the fact that the same correlation works for both 
young and millisecond pulsars and that it does not require 
enhancing the magnetic field using multi-polar components 
near the sur face of the star as wa s invoked in outer gap 
models (e.g., IZhang fc Cheng! (l2003l H are a new outcome of 
this study, even assuming physical content to the outer gap 
computation of /. 


4.4 Correlations with PLECl parameters 


In this section we study whether our best-fitting SC model 
parameters present any correlation with the PLECl phe¬ 
nomenological fits of the 2PC, as defined by 


-K F 
dE ^ 


-r2p 


exp 


— ) 

E/cutoff / 


( 8 ) 


Thus, here we focus on r 2 pc, Ecutoff, and the 0.1-100 GeV 
luminosity L, related to the normalization ifpc and to the 
distance (see Appendix|B]for details). We recall that the for¬ 
mer parameters, different from the set xq, Nq and Ey, come 
from a likelihood fit to the 7 -ray data, and are not associated 
to the predictions of any concrete theoretical model. 

To put any arising trend in persp ective it is interest ing 
first to use the PLECl fits from 2PC llAbdo et al.ll20l3) to 



















Fits of pulsar 'y-ray spectra 11 




2.0 3.0 4.0 


log(E,„„ff(MeV)) 



0.0 1.0 2.0 3.0 


r 


2pc 


Figure 5. Search for correlations between the three parameters of the best-fitting PLECl models, taken from the 2PC. Color coding is 
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see if there is any self-correlation among them for the pulsars 
in our sample. This is explicitly shown In Fig. [3 Whereas the 
luminosity distinguishes among YPs and MSPs, as we have 
commented above, neither P 2 pc nor the cutoff energy i?cutoff 
are able to clearly separate the sample. For instance, the low¬ 
est and highest values of F 2 pc are attained both by a YP and 
a MSP. The T-ifcutoff plane, particularly for MSPs, show an 
apparent positive trend, where the highest luminosities are 
achieved at the largest Fcutoff, what can be understood in 


the framework of a SC model as explained below when we 
discuss the plane L-E\\. 

4 . 4.1 Energy cutoff EcutoB 

The energy cutoff traces the radiation released by the most 
energetic particles, i.e., the ones that have reached the max¬ 
imum values of Lorentz factor. Thus, it is not surprising to 
find that the larger it is E^^ the larger are the Flcutoff values, 
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Table 4. Correlation fits between couples of well-constrained pa¬ 
rameters (solid circles in the figures). For the Alo-parameter all 
considered pulsar fits are well-constrained. We show the best lin¬ 
ear fits y = A + Bx, with the associated Pearson coefficients, 
r. Units of parameter values are the same used throughout the 
paper, and are omitted here. Only very significant correlations 
are reported (r > 0.85). Number in parenthesis give the pulsars 
considered. 


X 

y 

A 

B 

r 

(24/81) 

log Sji 

log XQ 

17.2±0.5 

-1.31±0.07 

-0.974 

log Sic 

log E|| 

4.95±0.08 

0.78±0.02 

0.916 

log Sic 

log XQ 

10.8±0.2 

-1.07±0.05 

-0.906 

(81) 

log Nq 

log L 

15±2 

0.64±0.06 

0.852 

YPs (14/59) 

log Ecut 

logs II 

-1.2±0.4 

0.59±0.06 

-0.881 

log £11 

log(a;o/Sic) 

5±1 

-0.9±0.1 

-0.921 

MSPs (10/22) 

log Eout 

logs II 

-1±1 

0.5±0.1 

-0.892 


for both sub-samples. What seems to be a rising relation is 
broken at the largest values i?||. It is intrinsic to the cur¬ 
vature losses to have a limit in the cutoff energies (see the 
saturated Lorentz factors in Fig[T]). No other clear correla¬ 
tions with Nq and xo are found. 

4-4-^ Power-law index V2^c 

The plane r 2 pc-£^|| (r 2 pc-a;o) might also present a hint for a 
possible positive (negative) trend (albeit not at a level pro¬ 
ducing r > 0.85). Both parameters and xo (the degen¬ 
eracy of many cases makes this obvious) play a role in reg¬ 
ulating the observed 7 -ray slope at sub-peak energies (and 
the spectrum in general). Finally, No does not present any 
obvious trend with the power-law index, which is expectable 
since there is also no relation between the 7 -ray luminosity 
and the slope (see below). 

4-4-3 y-ray luminosity L 

The luminosity in 7 -rays is a separating parameter (as we 
have seen above) between the two sub-samples of pulsars. 
Note that for some YPs, we only have upper limits due to the 
unknown distance (its upper limit is taken as the edge of the 
galaxy in a given direction). The 7 -ray luminosity relation 
with the normalization No is expected by definition. Hints of 
a correlation between L and if|| could be argued (but not at 
the level of r > 0.85 for the well-constrained pulsars in the 
sample). If further data confirms it, can be physically moti¬ 
vated by the fact that the largest Lorentz factors attained 
by particles In larger electric fields produce a more intense 
and more energetic photon flux. Since the radiation losses 
in MSPs are more efficient (closer light cylinder), particles 
are less energetic and the two sub-sample are divided in this 
plane. No clear trend is seen with xq, due to the larger errors 
in the determination of the latter parameter. 


4.5 Lightcurve parameters 

We have also tested the possible appearance of correlations 
among the three SC model parameters with lightcurve fea¬ 
tures, with the latter being as described in the 2PC. In par¬ 
ticular, we considered the number of peaks in the lightcurve, 
the radio lag, and the 7 -ray peak separation. In principle, 
one could expect some correlation of some of these param¬ 
eters with Xo or xo/Ric, if the latter values were effectively 
prompted by some geometrical effect. However, we find no 
significant trend in any of the planes. 


5 CONCLUSIONS 

This work represents the first systematic spectral fitting of 
7 -ray pulsar spectra using a physically-based model includ¬ 
ing particle dynamics along a magnetic field line and its as¬ 
sociated SC radiation. Our models are effectively described 
by only three parameters. In order to determine their value 
we have extensively explored thousands of models for each 
pulsar, calculating the evolution of the particle energetics 
and the locally emitted radiation in each. 

In a few cases on which we have individually commented 
above, residuals are seen above ~ 10 GeV. This can be 
due to a relatively important contribution of inverse Comp- 
ton scattering (e.g.^ as in the Crab and Vela pulsars, see 
Ahdo et al.ll2010'a H: iLvutikov et al.ir2012l : IVigano fc Torre^ 
2 OI 5 I for discussions), or been spuriously generated due to 
the superposition of quantitatively different phase-resolved 
spectra, all compatible with SC radiation (like in Geminga, 
IVigano fc TorredboiSh . 

For most of the pulsars studied, the best-fitting mod¬ 
els are satisfactorily reproducing the phase-averaged spectra 
in the whole Fermi-LAT range. These models are able to 
systematically explain the variety of both the energy peak 
and the flat slope of the spectra below ~ 1 GeV. The lat¬ 
ter property is often observed, and is not compatible with 
models where the spectra is generated only from the most 
energetic particles, i.e., when they have already lost their 
perpendicular momenta. Our results stand for the SG pro¬ 
cess as the main generator of 7 -ray pulsed emission, from 
most objects detected. We stress that pulsars can be fit by 
our models with only three free parameters, two of which 
are highly correlated. 

This study allow us to infer physically-motivated pa¬ 
rameters, i.e., the parallel electric field in the gap F||, the 
lengthscale over which most of the observed photons are pro¬ 
duced Xo, and the number of particles responsible for such 
detected emission No- Although our model does not include 
geometry, these parameters give an idea of the effective val¬ 
ues of these quantities in the region of the magnetosphere 
where the radiation is produced. Other parameters in the 
model have been neglected, being fixed to one and the same 
fiducial value for all pulsars, since they have a much mi¬ 
nor im pact on the trajectories and spectra (IVigand fc Torre^ 
I 2 OI 5 I) . This is the reason why we are not able to place any 
strong constraint on the location of the acceleration region, 
being able only to exclude the near-surface region. 

Some qualitative features are discovered by this system¬ 
atic approach. The accelerating electric field is typically of 
the order of 10® V/m for YPs, and 10® V/m for MSPs. There 
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is a strong anti-correlation between _E|| and xq, which is con¬ 
tinuous across pulsar types. The values of xo range from a 
few to hundreds of km, which represent a small fraction of 
the corresponding magnetospheric sizes (xo/i?ic ~ 10 “"*- 
10 “^). 

There is a strong correlation linking iJ|| with the mag¬ 
netic field at light cylinder, Bic, which is to be considered as 
a proxy to the magnetic field in the gap. This, together with 
the absence of any correlation with the surface magnetic 
field Bs, confirms that the place of high-energy emission is 
located in the outer magnetosphere of pulsars. Such B^f-Bic 
correlation quantitatively unifies under the same trend the 
two sub-classes, YPs and MSPs, which are otherwise (except 
in the plane i?||-xo, as commented above) seen as disjoint 
populations. In this paradigm, is larger in MSPs because 
their light cylinders are much closer to the surface than they 
are for YPs. 

In order to produce gamma-ray spectra up to GeV, 
MSPs need an accelerating electric field that looks large in 
comparison to their relatively low surface magnetization. In 
the usual OG m odel, and in order to g enerate such par¬ 
ticle acceleration, I Zhang fc Ghenel (120031 1 invoked multipo¬ 
lar magnetic fields existing close to the surface, about 1000 
times larger in magnitude than the dipolar surface value. 
Apart from the many assumptions that c aveat this gap 
model as a whole (see IVigano et al.l (l2015bll for a detailed 
discussion), the solution found here is simpler: corre¬ 

lates with the local value of B in the outer regions of the 
magnetosphere, i.e., close to or beyond the light cylinder. 
Whereas in previous studies, £'|| was inferred by adopting 
an over-simplihed estimation involving P and Ba, here i?|| 
is inferred by considering the physical process of electric ac¬ 
celeration, radiative losses, and the subsequently radiated 
photons. The assumptions (and these latter processes them¬ 
selves) make no qualitative difference between particles in 
the surrounding of MSPs and YPs. 

The systematic finding of low values of xq in comparison 
with Ru for both classes of pulsars may point to something 
fundamental related with the distribution of the particles. 
Although our modelling relies on an effective approach, the 
results indicate the need for a dominant number of particles 
with relatively low values of Lorentz factor, one or two orders 
of magnitude lower than the value of saturation P ~ 10 ^- 
10®. On one side, one could appeal to the cascading pro¬ 
cess as the mechanism able to replenish a great number of 
particles which quickly lose their perpendicular momentum. 
Second, the beaming effect, which scales with l/T favours 
the detection of low-energy photons. These effects, possi¬ 
bly responsible for the low-energy flatness of energy spectra 
in many pulsars, is, in our model, simplihed by adopting a 
larger weight to the initial part of the trajectories, i.e., small 
Xo - Numerical simulations involving cascading and geometry 
are needed to physically reinforce our interpretation, and to 
allow a comparison with pulse profiles. 

We also note that, both in radio and in 7 -rays, the gross 
spectral properties, like the luminosity, or the energy peak, 
show no or weak correlations with timing properties. Thus, 
although correlations between rotational energy and lumi¬ 
nosity have been proposed many times, we remind that the 
observations show a huge dispersion, which cannot be ne¬ 
glected. This strengthens the importance of the correlations 
we find, particularly because they unify YPs and MSPs. 


An important prediction is that, if we extrapolate the 
correlation Ejj(Bic) to the magnetar regime (Bic ~ 10 “^- 
10^ G due to their long periods, P ~ 2-10 s), it is clear that 
i?ll would be very weak, E^\ < 10 ® V/m, and unable to pro¬ 
vide particles energetic enough as to emit 7 -rays. Thus, we 
expect magnetars not to be visible cis pulsed 7 -ray emitters 
(unless a very different mechanism from SC radiation is at 
play). 

Finally, we note that with the forthcoming publication 
of the third Fermi-LAT pulsar catalog, with a longer time 
baseline and the new Pciss 8 analysis, a larger and better 
quality sample could be analyzed. With such sample, for 
which more than 1/3 will consist of pulsars not contained in 
our study, we shall look deeper and in an independent way 
into the trends here illustrated and their physical implica¬ 
tions. 
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APPENDIX A: DETAILS OF THE FITTING 
MODEL 

A1 Formulae of SC radiation 

The sinKle-particle SC ra diation power is dCheng fc Zhan j 
ll99fil : IVigano et alJl2015al l 

= (Al) 

where 


z = (Q2reff) ^ , 

POO 

P{y) = / Kz/3{y)dy', 

J y 


V = 


_E 

'K ’ 


E, = l&Qzr®, 

m,c?V sin a 


eB 


Q 2 — 


cos^ a 


1 + 3C + e + 




rc sm a 


Tgvr cos^ a 


refi 


cos^ a 
2 


l + ^ + 


gyr 

Tc 

-2 


^ rj [l + 7(reffQ2)-"] 

rig 8(Q2refi)-i 


(A2) 

(A3) 

(A4) 

(A5) 

(A6) 

(A7) 

(A8) 

(A9) 

(AlO) 


Here, m and F are the rest mass and the Lorentz factor of 
the particle, a, and are the pitch angle, the Larmor 
radius, and the radius of curvature of its trajectory, respec¬ 
tively, e the elementary charge, B the local strength of the 
magnetic field, h the reduced Planck constant, c is the speed 
of light, Kn are the modified Bessel functions of the second 
kind of index n, E is the photon energy, Ec is the charac¬ 
teristic energy of the emitted radiation. In the limits of high 
or vanishing perpendicular momentum, Eq. reduces to 
purely synchrotron 3> 1) or curvature radiation <C 1) 
formulae, respectively. In any case, the peak of the spectrum 
is located close to Ec and, for energies E E^, the domi¬ 
nant term in Eq. (lAlll . F{y), provides dPsc/dE ~ E^'^^. 


A2 Particle dynamics 


For each model, we simulate the motion of charged parti¬ 
cles by evolving the parallel and perpendicula r momenta of 
partic les according to the equations of motion dVigano et al.l 
l2015all 


d{p sin a) 
dt 

d{p cos a) 
dt 


Pac sin a 
V 

Psc cos a 
eE||- 

V 


(All) 

(A12) 


where p = Fmt; is the momentum of the particle, v its spa¬ 
tial velocity, E|| is the parallel component of the electric 
field, and Psc is the single-particle SC power, obtained by 
integrating Eq. (EH in energy: 


P — 

sc - 


2e^r^c 

3r2 “ 


(A13) 


At each position in the gap we use the corresponding lo¬ 
cal values of B{x) = Ba{R*/x)*‘ and rc{x) = Ric{x/Ric)^ 
(b and rj ar e model parameters, here fixed to 2.5 and 0.5, 
respectively, IVigano et al.ll2015a : IVigano fc Torresir2015l ~l , to 
compute Eqs. (IAlll) - (IA12ll consistently. 

In Fig. [T] and the associated text, we describe the main 
outcome of our simulated trajectories. Here we derive why 
the correlation Ey ~ l/io, discussed in the text and shown 
in many cases in Figs. IB6IB101 is expected to appear. When 
the electric acceleration term dominates over the SC power, 
eE|| Psc cosa/v, one can simplifies Eq. (IA12II . considering 
also that cos a —>■ 1 very fast, as seen in t he numerically com¬ 
puted trajectories dVigano et al.|[^15al l: dF/dt ^ eE^^fmv. 
Since t = x/v, then, in the accelerating regime, Ey ~ I/®. 
In other words, the linear rise of the trajectory (r(a:), left 
panels of Fig. [TJ can be displaced to lower x by rising the 
value of E||. In many pulsars with relatively flat slopes, this 
part of the trajectory is dominating the spectrum, and the 
best-fitting models are endowed with a value of xo smaller 
than the distance at which saturated values of F are reached. 
Therefore, the best-fitting solutions of these particular pul¬ 
sars will have an intrinsic degeneracy: larger Ey can be com¬ 
pensated by re-scaling xo with 1/Ey. This quantitatively ex¬ 
plains the correlation. 

As discussed in depth in tj4.1[ such analytical argument 
explains the individual degeneracy appearing in some of the 
Figs. IB6IB101 but not the similar global trend found when 
comparing the best-fitting parameters for all pulsars, i.e., 
again, Ey ~ 1/xo (top panels of Fig.[2]and first row of Table 

H). 


APPENDIX B: FITS OF THE SAMPLE 

Figures of this Appendix show individual best-fitting models 
compared with pulsar data, and contours of x^/Xmin in the 
plane logEy — \og{xo/Ric)- In each of the contour plots, 
we mark with a white cross the location of the best-fitting 
model. The extreme values of xo, Ey and No within the 
region having IXmiu < 1-3 are taken as the errors of model 
parameters and are quoted in Tables [5] and [3] (see also l3.2.51) . 
Due to the finite steps of the explored parameter values, 
the resolution of the contours can be visually insufficient at 
times, especially when data suffer very small relative errors 
(like in the Vela pulsar). Red and blue lines or labels indicate 
YPs and MSPs, respectively, following the color coding used 
in the main text. 

Observational spectra are taken from the publicly avail¬ 
able data of the 2PC. Data consist of a number of Abin ~ 
5—10 energy bins (with extremes Ei and E 2 , and weighted 
central energy Ecent), each one with its associated photon 
flux, E^JOs (in units photons cm“^s“^), and its associated 
statistical error, (5E,(b". We consider the isotropic luminos¬ 
ity, 

Aoi," = (Bl) 

where d is the distance of the pulsar to the Earth, and we 
neglect possible beaming effects which would reduce the in¬ 
ferred luminosity. We plot the binned functions E^dN/dE = 
A'centT/'^'", in units erg s“^, both for the theoretical models 
and data. 

To explore the space of parameters, we span a grid of 
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101x51 equi-spaced values of [log-By, log xo], covering two 
and five orders of magnitude, respectively. The range of By 
varies case by case, while xo always spans [0.0001 — lOjBiC. 
We also evaluate the case of uniform effective distribution 
of particle. For each pair of values (By, so), we evaluate the 
expected spectrum over a grid of hundreds of points in the 
0.1-100 GeV range. Then, we renormalize it by the best- 
htting value of Nq , found by means of scanning a grid having 
progressively finer steps, up to the numerical convergence, 
defined by SNo/Nq < 10“®. Therefore, we integrate the SC 
photon spectrum (related to the SC energy spectrum, Eq. [Jl 
by dA^gap/dB = (l/B)dPgap/dB) in each bin (and normal¬ 
ized by the bin width), to obtain the binned SC photon 
spectrum (number of photons per unit energy): 



(B2) 


With this we calculate the goodness-of-fit indicator: 



(B3) 


where dBobs = 47rd^(5B,5i(“ is the luminosity error for each 
bin. Such error neglects the (likely large) uncertainties on 
distance and beaming factor, but this would only change 
the luminosity, thus providing a different (but accordingly) 
inferred normalization Nq, i.e., without changing the spec¬ 
tral shape, and, leaving the constrains on By and xo / Bic -as 
well as all trends- unscathed. When distance is unknown, we 
use its upper limit (the edge of the galaxy in that particular 
direction), and the related upper limit on B, to constrain Nq 
(which, therefore, has to be taken as an upper limit). 
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Figure Bl. Fits of the Fermz-LAT spectral data and models for the pulsars considered in the sample (I). We plot the binned functions 
E^dN/dE = , in units erg s“^, both for the theoretical models and data. Red and blue lines indicate YPs and MSPs, respectively. 
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Figure B2. Fits of the Fermz-LAT spectral data and models for the pulsars considered in the sample (II). 
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Figure B3. Fits of the Fermz-LAT spectral data and models for the pulsars considered in the sample (III). 
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Figure B4. Fits of the Fermi-LAT spectral data and models for the pulsars considered in the sample (IV). 
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Figure B5. Fits of the Fermz-LAT spectral data and models for the pulsars considered in the sample (V). 
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Figure B6. X^/Xmin contours on the log£^|j-log(a:o//?ic) plane for the pulsars considered in the sample (I). Red and blue labels indicate 
YPs and MSPs, respectively. 
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Figure B7. X^/Xmin contours on the log£’|j-log(a:o/-Ric) plane for the pulsars considered in the sample (II). 
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Figure B8. X^/Xmin on the log£^|j-log(xo//lic) plane for the pulsars considered in the sample (III). 
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Figure B9. X^/Xmin contours on the logF||—log(a:o/i?ic) plane for the pulsars considered in the sample (IV). 
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Figure BIO. X^/Xmin contours on the logF'ij—log(xo/-Ric) plane for the pulsars considered in the sample (V). 
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